Many-body Correlation Effect on Mesoscopic Charge Relaxation 
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We investigate in a nonperturbative way the dynamics of a correlated quantum capacitor. We find that the 
many-body correlations do not disturb the universal low-frequency relaxation resistance per channel, Rq (iu = 
0) — h/4e^ ensured by the Korringa-Shiba rule whereas the interpretation of the quantum capacitance Cg 
in terms of the density of states fails when strong correlations are present. The AC resistance Rq{ijj) shows 
huge peaks (with values larger than h/Ae^) at fitj f» ±r*, where F* is the renormalized level broadening. 
These peaks are merged to a single one at a; = when a finite Zeeman field is applied comparable to V* . The 
observed features of Rq, being most evident in the Kondo regime, are attributed to the generation of particle-hole 
excitations in the contacts accomplished by spin-flip processes in the dot. 
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PACS numbers: 73.63.-b, 74.50.+r, 72.15.Qm, 73.63.Kv 

Introduction. — Optoelectronic devices such as Ught 
sources emitting single photons on demand are of an enor- 
mous interest in quantum information. Recently, its solid-state 
analogue, a quantum capacitor (QC), was created by Feve et 
at. Ill], in which a quantum dot (QD) was coupled to a single 
reservoir via a quantum point contact. Fast time-controlled 
variations of fractions of nanoseconds on the dot gate poten- 
tial produced the single-electron source of emitting electrons 
in a desirable quantum state. Previously, Gabelli et al. [2] 
showed that the QC could act as a RC circuit with a quan- 
tized resistance as predicted in a series of seminal works [3,- 
[st]- The experimental interest in the AC properties of meso- 
scopic conductors |6-11] has been revived recently due to 
the experimental confirmation of the quantization of both the 
AC current HI] and the quantum resistance ]2]. In a coherent 
conductor the AC transport is highly sensitive to the internal 
distribution of charges and potentials that need to be calcu- 
lated in a self-consistently manner to ensure a gauge invariant 
and cuiTent conserving description. For a macroscopic capaci- 
tor the low-frequency dynamical conductance depends on two 
elements: the geometrical capacitance C, and the resistance 
R. Classically, the electric field on the surface of the metallic 
plates is completely screened and C is characterized solely via 
Coulomb forces. However, Biittiker et al. pointed out that for 
a coherent nanoscale system that rule statement is not longer 
valid, since electric fields penetrate at distances of the order 
of the Thomas-Fermi screening length ^ 0, [12. [HI] ■ Con- 
sequently, the capacitance for a coherent conductor, termed 
as electrochemical capacitance, depends on the geometry, on 
its physical properties and particularly on its density of states 
(DOS) through the quantum capacitance Cg. Even more sur- 
prisingly, the resistance becomes quantized independently of 
the value of the transmission through the mesoscopic con- 
ductor 

To distinguish it from the DC resistance, 
Rq is called charge relaxation resistance, and together with 
= (C-i + C'-^y^ defines the RC time for a QC, that is. 



the charge relaxation time upon the action of time-dependent 
potential. 

So far several aspects of AC conductance for RC quantum 
circuits have been addressed 11 1 8142511 focusing mostly on the 
spin polarized case. Therefore the understanding on the ef- 
fect of many-body correlations in RC circuits is still missing. 
This Letter attempts to fill this gap, offering a physical pic- 
ture of the influence of many-body interactions on the dynam- 
ics of QCs. We consider a QC formed by an ultrasmall QD 
and calculate its Cq and Rq by using their relations with the 
charge susceptibility in the linear response regime. We per- 
form a thorough study by tuning the QD from the resonant 
tunneling regime to the Kondo regime. In the Kondo regime 
virtual tunneling transitions between the dot and the reservoir 
flip efficiently the dot spin resulting in a many-body singlet 
spin state with binding energy given by the Kondo energy 
UbTk- First, we find that the Korringa-Shiba (KS) relation 
02611 . valid only in the Fermi-liquid regime, ensures the quan- 
tization of Rq even in the Kondo regime as long as the fre- 
quency is low enough, hbj ^ UbTk. However, it is pointed 
out that the relation of Cq with the localized DOS should be 
revised when the charge dynamics is frozen due to Kondo cor- 
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FIG. 1: (color online) (LEFT) Illustration of a quantum capacitor 
and (RIGHT) second-order tunneling processes that generate a single 
particle-hole pair in the conduction band without [(a)] and with [(b)] 
a spin flip in the dot. Filled and empty arrows in the contact indicates 
particles and holes, respectively. Here we assume > 0. 
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relations. More importantly, the AC resistance Rg (ui) displays 
peaks at huj ^ ±r*, not present in a mean-field description 
for the Coulomb interaction. Here F* is the renormalized level 
broadening of the major tunneling dot level and F* — kgTji 
in the Kondo regime. Third, a departure from the universal 
value for Rq{u) ~ 0) is achieved when an applied Zeeman 
field IS.Z is comparable to F*. The enhancement of Rq at zero 
and finite frequencies is attributed to the creation of particle- 
hole ip-h) excitations accompanying a spin-flip in the dot due 
to the strong Coulomb interaction, see Fig.lifb). 

Model. — The essential features of a mesoscopic RC cir- 
cuit in the presence of correlations can be well captured in the 
Anderson model where an interacting single-level QD is cou- 
pled to a single-channel electron reservoir subject to a weak 
time-dependent voltage V{t) = V^ccoswt. The Hamiltonian 
of this system is H = Hl + + ^t- The lead part, 
— Z^k/ji^k + eV^(t)]cj^^Ck^, describes the noninteract- 
ing conduction electrons with energy et (measured with re- 
spect to the Fermi energy ep ^ 0) and spin fi in the reservoir, 
and the tunneling of electrons between the reservoir and the 
dot is modeled by Ht — X^k^i [^kdj^Ck^ + (/i.e.)] in terms 
of energy-independent tunneling matrix element, tk ~ t. 
The hybridization between the dot and the lead is character- 
ized by a tunneling amplitude F = 7rpo|ip iPo is the con- 
tact DOS at the Fermi energy). The dot Hamiltonian reads 



-2i?cf^t"i' where = dj^d^ is 



the dot occupation operator and Eq — e /2C is the Coulomb 
charging energy. The orbital level = — ij,Az/2 is spin- 
dependent due to a Zeeman energy Az- The time-dependent 
voltage V{t) induces the polarization charges Nij{t) between 
the dot and the gate, which in turn leads to the time-dependent 
potential C/(t) — |e|A^;7(<)/C inside the dot. Consequently, 
the applied voltage not only generates a current I{t) between 
the lead and the dot, but also induces a dot-gate displace- 
ment current Id{t) = e{dNu/dt) = -C{dU/dt). Charge 
conservation requires /(<) + Id{t) — 0. Assuming that the 
gate-invariant perturbation V{t) — U{t) is sufficiently small, 
the linear response theory leads to the relation, /(w) — 
g{oj){V{oj) - U{uj)), where g{t) = (ie/h) {[I{t),J\f])e{t) 
is the equilibrium correlation function between the occupa- 
tion operator M — '^^'^ '^^e current operator I = 
e{dAf/dt). Note that the current-density correlation function 
g{uj) is directly related to the charge susceptibility Xc{t) — 
—i {[N'{t),Af]) &{t), which is preferable for numerical com- 
putation, via the relation g{uj) = iui{e^ /h)xc{^)- Then, with 
the help of I{uj) = —Id{i^) = —iLuCU{uj), the dot-lead 
impedance Z{u!) — V{uj)/I{uj), which is experimentally ac- 
cessible, is given by Z{uj) = l/(— iojC) + l/g{uj). The re- 
laxation resistance and the quantum capacitance are then ex- 
pressed in terms of the charge susceptibility as 



= Re 



27riwxc(w) 



Cqiuj) 



= Im 



27riXc(w) 



(1) 



The numerical normalization group (NRG) 027112811 treats the 
Coulomb interaction in a nonperturbative way being the most 
adequate method for computing the charge susceptibility L29l1 . 
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FIG. 2: (color online) (a) Zero-frequency limits of Rq and Cq versus 
Ed. (b) Typical spectral structure of the real and imaginary parts of the 
charge susceptibility Xc (c,d) Enlarged views of Im[xc] near w = 
for two values of e^- The dashed (blue) and dotted (green) lines are 
tangent lines at infinitesimally small a; and at points (indicated by 
arrow) where the slope is maximal, respectively, (e, f) Rq (to) in the 
logarithmic scale. Here F = 0.04, and Ec = 0.5. 



While the imaginary part of the susceptibility is obtained by 
the NRG procedure its real part is calculated via the Kramers- 
Kronig relation. Note that the NRG results for the finite- 
frequency linear response in the Kondo regime are known to 
be reliable as long as the perturbation is weak enough llsoll . 
We focus on the zero-temperature case and use the contact 
bandwidth D as the energy unit. We set fcs = 1 hereafter 

No Zeeman splitting, Az — 0. — Figure |2 a) shows our 
main results for Rq and Cq for the spin-degenerate case. 
First, the zero-frequency limit of the relaxation resistance, 
Rq{Ld — )■ 0) is always close to the universal value h/Ae^, re- 
gardless of values of and Ec- This value can be interpreted 
as the composite resistance of two parallel resistors of resis- 
tance /i/2e^, the well-known universal resistance per channel 
llsl-Isi l20[l ■ The NRG results show the quantization of charge 
relaxation even in the Kondo regime where many-body corre- 
lations are effective. The observed small deviations from the 
exact value h/Ae^, persisting even in the noninteracting case, 
are attributed to the finiteness of the contact bandwidth D, 
which introduces a frequency-dependent real part into the dot 
self energy , Re[i;(a;)] = ~{T /2TT)\n\{D - hw) / {D + nijj)\. 
Its presence slightly violates the Fermi-liquid assumptions 
and, consequently, the KS relation is not exactly fulfilled so 
that the universal value is not recovered: Since Re[I](ti;)] in- 
creases with uj in magnitude, the deviations are larger as the 
resonant level becomes far from = 0. The universal value 
can be restored by setting all the relevant energy scales to be 
much smaller than D [31]. Second, the quantum capacitance, 
Cq exhibits two remarkable considerations: (i) at the degen- 
erate points, €d ~ cf and + 2Ec ^ ep, Cq shows two 
pronounced peaks [see Fig. |2a)], which is consistent with 
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the known understanding that Cq is proportional to the dot 
DOS Pdoti^p) iSBIl^l. and (ii) Cq remains quite small in 
the Kondo regime although the DOS at the Kondo resonant 
level pinned at the Fermi level achieves its maximum value. 
It implies that the Kondo resonant level, even though it can 
open a tunneling channel, is not a real level which can hold 
real charges and cannot contribute to the capacitance. Hence, 
in the presence of many-body correlations Cq is not, always, 
directly related to the DOS. 

The frequency dependence of Rq{uj) and Cq{uj) is analyzed 
in Fig. |2] where the the behavior of the real and imaginary part 
of the charge susceptibility Xc{^) is shown, see Fig.|2b) for 
Cd = —0.3. The imaginary part of Xc(w) reflects the cou- 
pling between the ground state and p-h excitations due to the 
dot-lead hybridization. Since the spectral density of multiple 
p-h excitations increases with energy, |Im[xc]| would grow 
monotonically with |w|. However, a finite D puts an upper 
limit to the energy for p-h excitations [\uj\ > 0{D)] resulting 
in the observed nonmonotonic behavior for Im[xc]- More- 
over, Im[xc] has two kinks at |aj| = min(|e(j|, \ed + 2£'c|) 
since beyond this frequency excitations accompanied with 
a charge excitation contributes to Im[xc] as well. An interest- 
ing structure appears in Im[xc] near a; = 0, see Figs. |2c) 
and (d) for two dot level positions, — —0.1 (correspond- 
ing to the fluctuating valence regime) and = —0.3 (Kondo 
regime). Close to a; = 0, Im[xc] depends linearly with lu, 
mainly due to single p-h excitations [see tangent (dashed) 
lines in Figs. |2|c) and (d)]. However, Im[xc] departs from 
linearity when ut becomes of the order of the effective hy- 
bridization r* (= Tk in the Kondo regime). The renormal- 
ized hybridization T*(Tk) is extracted from the width of the 
resonance close to(at) ep in pdot- Besides, we found that the 
slope of Im[xc(w)]/w is the largest at hu ^ F*, while the 
change in Re[xc(w)] is marginal. As a consequence, Rqicd) 
(see Eq. ([T]i) exhibits two side peaks at huj = ±r* as shown 
in Figs. |2e) and (f). Notice that in the Kondo regime Rqicu) 
becomes much larger in order of magnitude than the univer- 
sal value, see Fig. |2jf). Remarkably, such peak structure in 
Rq{oj) is absent in the noninteracting case. For a noninteract- 
ing system, the analytical expression of Rq{uj) in the wide- 
band limit is given by L32J i?r"(w) = {h/e^){[G{uj) + 
(fta;/r)F(cj)]/[G(cj)2 + F{u}y]} wifli G{uj) = ln{[(e<j + 

+ r2][(ed - f^? + + r']} and Fiio) = 

2{tan-i[r/(ed - hio)] - tan-i[r/(erf + huj)]}. 
increases monotonically with increasing \uj\, and the only 
characteristic energy scale is Cd- Hence, the peaks seen in 
Figs. |2e) and (f) are a genuine many-body effect. In sum- 
mary many-body correlations, apparently having no impact 
on the zero-frequency value do affect Rq{uj) at finite frequen- 
cies, by forming a pronounced peak at hui — r*(Tx). The 
explanation for these observations will be given later 

Finite Zeeman splitting, ^ 0. — The spin-split case 
(A^ > 0) in the presence of external magnetic fields is il- 
lustrated in Fig.|3] Interestingly, — 0) versus A^ ex- 
hibits a peak structure reaching values much larger than the 
quantized resistance in the spin-degenerate case: for example. 
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FIG. 3: (color online) (a,b) Rq{u> — 0) and Cq{uj — 0) versus Az 
for annotated values of ed with Ec ~ 0.2. The enlarged view of 
Rq{ijj — 0) for large Az is shown in the inset, (c) Contour plot of 
Rq as a function of cj and Az in the logarithmic scale in the Kondo 
regime (e^ = —0.1). (d) Rq{io — 0) versus Az for different values 
of Ec as annotated. Here F = 0.02 is used. 



Rq{u! = 0)|max ~ 100 X /i/4e^ for = -0.15. Furthermore, 
the peak is exactly located at A^ = T*{Tk) for the fluc- 
tuating valence (Kondo) regime. The peak height increases 
as the the effective hybridization decreases so it is the high- 
est in the Kondo regime. In the meanwhile, Cq{uj 0) re- 
mains rather constant, except at the resonant tunneling regime 
(e^ w 0) where it displays a small peak, see Fig. [3jb). The 
evolution of the spectral distribution of with Az is dis- 

played in Fig.[3lc) for — —0.3 (Kondo regime). As Az 
increases, the low-frequency part of Rq{ijj) for \hiLi\ < ksTx 
keeps going up until Az reaches ksTK', the side peaks are 
merged into the central peak. With increasing Az further, 
the central peak diminishes gradually and, eventually, together 
with the side peaks located athw — ztknTK, disappear com- 
pletely. We have observed a similar transition of Rqioj) with 
Az in the resonant tunneling regime (e^ = —0.05, 0) ex- 
cept that the variation of the central part is smaller. Finally, 
Fig. [3jd) compares the zero-frequency values of Rq for dif- 
ferent values of the Coulomb interaction in the resonant tun- 
neling regime. In the noninteracting case, there is no peak at 
all, with Rq{uj — s- 0) equal to h/Ae^. However, as soon as the 
charging energy 2Ec becomes comparable to F* ^ F, a peak 
starts to rise up and manifests itself for 2Ec 3> F. It implies 
that the existence of the peak structure observed in Fig. [3|d) 
definitely has its origin in the Coulomb interaction. 

Discussion. — Now we have two questions to be an- 
swered: (l) How can Coulomb interaction increase the re- 
laxation resistance far beyond the universal value, h/Ae^ and 
(2) Why does it take place noticeably at huj = ±r*{TK) 
in the fluctuating valence (Kondo) regime for Az = or at 
hoj = for Az = F*(rR-)? The charge relaxation resistance 
is attributed to p-h pair generation in the conduction band as 
shown in Fig. [T] Such processes are put in action when the 
dot-lead tunneling is switched on. The tunneling in turn hy- 
bridizes dot and conduction band electrons, resulting in low- 
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ering of the ground state by the effective binding energy T* 
(Tk in the Kondo regime). It means that the p-h generation 
starts when the energy suppHed by the source is larger than 
r*. This argument explains the observed peak in Rq{oj) at 
hoj — ±r* in the absence of the Zeeman splitting. In the 
presence of finite but small Zeeman splitting, the energy cost 
can be compensated by the Zeeman splitting. The p-h pair 
excitation states shown in Fig.fljb) are now lowered by Az, 
and when Az « T* they become almost degenerate with the 
ground state, allowing p-h pair generation with negligible en- 
ergy cost. Hence, Rqicd) exhibits a single peak at hu = 
when Az — F*. This argument works solely when Az < T* 
in which the ground state is not yet completely polarized and 
there exists a finite coupling among spin-down dot states and 
spin-up dot states accompanying with a p-h pair generation in 
the reservoirs, see Fig.[T] The importance of the spin flip in the 
boosting of the relaxation resistance also explains why Rq{oj) 
can reach higher values in the Kondo regime. The Kondo 
ground state is built from spin fluctuations due to spin-flip 
scattering among the localized dot electron and the delocal- 
ized electrons in the reservoirs, thus spin-flip processes have 
large amplitudes in its wavefunction. Hence the processes as 
shown in Fig. [Tib) can happen more frequently, leading to a 
large Rq. Similarly, the answer for the first question is now 
ready. The spectral weight for the charge correlation function 
is proportional to | (ajA/'Igs) p, where a represent the excited 
states. In the second-order perturbation theory, this weight 
corresponding to the processes in Fig. [Tib) is given by 



{a\M\gs} 



^ Ay. \ E„ e„ 



(2) 



in the w — ?► limit with = 2i?c + e^i - Interestingly, this 
weight vanishes for Eq — for any value of Az- Thus, for 
the noninteracting case there exists no p-h pair generation pro- 
cess accompanying spin flip in the dot, and no boosting of the 
relaxation resistance can happen. For finite values of Eq, the 
weight is finite [see Eq. (O] and for Eq, oo, it becomes 
t'^ I [t^t^)"^ . This value can be substantial depending on the 
level position. Note that this analysis is not correct quanti- 
tatively because high-order events should be considerably in- 
volved in the observed phenomena. The observed boosting of 
Rq at Az ^ r* indicates that the perturbation in the dot-lead 
tunneling or F is risky. A more general theoretical analysis 
that treats Az and F on equal footing could provide more 
quantitatively reliable interpretation. Besides, this perturba- 
tive analysis does not work in the Kondo regime where the 
strong dot-lead coupling is important. One may want to study 
the Kondo regime by an effective single-particle Hamiltonian 
with a dot level at the Fermi energy with the effective hy- 
bridization Tk- However, this picture is only suitable in the 
Fermi-liquid regime in which p-h excitations accomplished 
by spin-flip events in the dot are not allowed. Besides, this 
effective model predicts an enhanced mesoscopic capacitance 
Cqiyj — > 0) due to the presence of the resonant level at the 
Fermi level. As noted before, however, the Kondo resonant 



level cannot contribute to the charging of real charges. 

Conclusion. — In closing, we have investigated the dy- 
namics of a many-body quantum capacitor. Using the rela- 
tion of charge relaxation resistance and quantum capacitance 
with the charge susceptibility, we find that in the deep Kondo 
regime the KS rule ensures a quantized Rq(uj = 0) = h/Ae^- 
Besides we show that the interpretation of Cq in terms of 
the DOS fails when many-body effects are present. Here, 
Cq{ijj = 0) becomes very small even when pdoti^p) becomes 
large due to a completely frozen charge dynamics in the pres- 
ence of Kondo correlations. Finally, we find that -Rg(w) is 
built by dot-lead tunneling events connecting p-h excitations 
in the reservoirs with spin-flip processes in the dot . This in- 
terpretation explains our results for Rq{Az = 0, uj) showing 
peaks at ftw = ±T*{Tk) and for R{Az = r*{TK),uj) ex- 
hibiting a peak at hoj ~ 0- 
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